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Stabilization  of  Rigid  Body  Dynamics 
by  Internal  and  External  Torques 

A.M.  Bloch*,  P.S.  Krishnaprasad**, 

J.E.  Marsden***,  and  G.  Sanchez  de  Alvarez**** 


Abstract 

In  this  paper  we  discuss  the  stabilization  of  the  rigid  body  dynamics  by  external  torques  (gas  jets)  and 
internal  torques  (momentum  wheels).  We  compare  the  stabilizing  quadratic  feedback  law  for  a  single  external 
torque  recently  analyzed  in  Bloch  and  Marsden  [1989b,c]  with  quadratic  feedback  torques  for  internal  rotors. 
We  show  that  with  such  torques,  the  equations  for  the  rigid  body  with  momentum  wheels  are  Hamiltonian  with 
respect  to  a  Lie-Poisson  bracket  structure.  Further,  these  equations  are  shown  to  generalize  the  dual-spin 
equations  analyzed  by  Krishnaprasad  [1985]  and  Sanchez  de  Alvarez  [1986].  We  establish  stabilization  with 
a  single  rotor  by  using  the  energy-Casimir  method.  We  also  show  how  to  realize  the  external  torque  feedback 
liquations  using  internal  torques.  Finally,  extending  some  work  of  Montgomery  [1990],  we  derive  a  formula 
for  the  attitude  drift  for  the  rigid  body-rotor  system  when  it  is  perturbed  away  from  a  stable  equilibrium  and  we 
indicate  how  to  compensate  for  this. 
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Introduction 

The  problem  of  stabilization  of  the  rigid  body  and  systems  of  rigid  bodies  is  of  importance 
for  numerous  practical  applications.  There  has  been  much  work  recently  on  stabilizing  the  angular 
momentum  equations  and  attitude  equations  of  the  rigid  body  with  n  >  2  torques.  Work  in  this 
area  includes  that  of  Baillieul  [1981],  Bonnard  [1986],  Brockett  [1976],  [1983],  Crouch  [1986], 
Aeyels  [1985a,b],  Krishnaprasad  [1985],  Sinchez  de  Alvarez  [1986],  [1989],  Aeyels  and 
Szafranski  [1988],  Sontag  and  Sussman  [1988],  and  Byrnes  and  Isidori  [1989]. 

A  related  area  where  there  has  also  been  much  progress,  is  the  problem  of  analyzing  the 
stability  of  coupled  rigid  and  flexible  bodies.  A  method  of  analysis  based  on  ideas  of  geometric 
mechanics  for  these  systems  (including  dissipation)  was  introduced  in  Krishnaprasad  [1985], 
Krishnaprasad  and  Marsden  [1987],  and  Baillieul  and  Levi  [1987].  In  particular,  stability  of  a 
rigid  body  with  flexible  attachment  was  analyzed  in  the  body  frame  by  the  Energy-Casimir  method 
developed  by  Arnold  [1966]  and  further  developed  by  Holm,  Marsden,  Ratiu  and  Weinstein 
[1985].  (The  method  has  a  long  history,  special  cases  of  which  can  be  found  already  in  the  last 
century  in  the  work  of  Riemann  and  Poincare  on  gravitating  fluid  masses).  In  this  method,  energy 
and  momentum  are  used  together  to  prove  Lyapunov  stability.  More  recently,  a  variant  of  this 
method  called  the  energy-momentum  method,  where  the  stability  analysis  is  in  the  material 
representation,  has  been  developed.  See,  for  example,  Simo,  Posbergh  and  Marsden  [1989]  and 
Simo,  Lewis  and  Marsden  [1989], 

Bloch  and  Marsden  [1989b,c]  showed  that  the  Energy-Casimir  method  could  be  used  to 
prove  a  stabilization  result,  namely  that  the  angular  momentum  equations  of  the  rigid  body  can  be 
stabilized  about  the  intermediate  axis  of  inertia  by  a  single  external  torque  applied  about  the  major 
or  minor  axis.  Moreover,  they  showed  that  there  was  an  interesting  Lie-Poisson  structure 
associated  with  the  system.  This  gave  a  quite  different  feedback  law  and  method  of  analysis  for  a 
closely  related  result  originally  due  to  Aeyels  [1985]. 

A  striking  feature  of  the  analysis  is  that  there  are  still  conserved  quantities,  even  while 
torque  is  being  applied.  This  leads  one  to  conjecture  that  the  feedback  might  be  realizable  as  an 
internal  torque  between  the  rigid  body  and  an  attached  rigid  body  or  bodies.  We  show  in  this 
paper  that  the  externally  stabilized  feedback  system  may  be  realized  as  a  rigid  body  with  three 
internal  rotors.  Our  analysis  is  based  on  the  original  analysis  of  this  system  by  Krishnaprasad 
[1985]  (see  also  Sanchez  de  Alvarez  [1986]).  In  fact,  we  show  that  the  three  internal  rotors  can 
realize  any  external  torque  feedback  for  the  rigid  body.  Moreover,  we  show  that  a  particular  choice 
of  internal  torques  for  body-rotor  system  makes  this  system  behave  precisely  like  the  classical 
heavy  rigid  body. 
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We  go  on,  however,  to  do  more  than  this.  We  analyze  the  rigid  body  with  internal  rotors 
with  certain  quadratic  feedbacks  that  make  the  rigid  body-rotor  system  Hamiltonian  with  respect  to 
a  Lie-Poisson  bracket.  We  show  that  with  a  specific  choice  of  feedback,  the  rotors  rotate  at 
constant  angular  velocity —  i.e.,  we  obtain  as  a  special  case,  the  driven  dual  spin  satellite  as 
analyzed  by  Krishnaprasad  [1985]  and  Sanchez  de  Alvarez  [1986]. 

Thus,  we  have  produced  a  further  class  of  systems  that,  despite  having  feedback  torques, 
are  Hamiltonian.  In  fact,  we  show  that  under  certain  integrability  conditions,  even  under  cubic 
feedback,  we  get  Hamiltonian  systems. 

We  then  discuss  the  Hamiltonian  structure  of  the  single  rotor  case  and  show  that  we  can 
again  prove  stabilization  of  the  system  about  the  intermediate  axis  for  sufficiently  large  torque  by 
using  the  Energy-Casimir  method. 

We  then  show,  using  the  work  of  Marsden,  Montomery  and  Ratiu  [1990],  Krishnaprasad 
[1990]  and  Montgomery  [1990],  (see  also  Levi  [1990])  on  geometric  phases,  how  an  attitude  drift 
can  occur  if  the  body  is  perturbed  away  from  a  stable  equilibrium.  Generalizing  a  formula  of 
Montgomeiy  [1990],  we  show  how  to  calculate  this  drift  precisely  and  indicate  how  to  compensate 
for  it. 

The  outline  of  the  paper  is  as  follows.  First  we  present  the  rigid  body  with  a  single  external 
torque.  We  introduce  feedback,  discuss  the  Lie-Poisson  structure  of  the  system,  and  recall  the 
stability  result  proved  in  Bloch  and  Marsden  [1989b,  c].  We  also  show  how  this  feedback  system 
is  indeed  the  generalized  (possibly  indefinite)  rigid  body.  Next,  we  show  how  to  realize  the  rigid 
body  with  external  torque  as  a  system  with  three  internal  rotors.  We  then  discuss  the  body-rotor 
system  with  quadratic  feedback  and  its  Hamiltonian  structure,  as  well  as  the  cubic  feedback  case. 
In  the  next  section  we  discuss  stabilization  by  one  rotor  with  quadratic  feedback  using  the  Energy- 
Casimir  method. 

Finally,  we  discuss  the  question  of  attitude  drift  and  how  to  compensate  for  it  We  remark 
that  the  work  discussed  above  on  phases  may  be  useful  here  when  used  in  conjunction  with  the 
work  on  chaos  in  mechanical  system  using  Melnikov's  method,  as  was  done  in  Oh  et.  al.  [1989]. 

Another  thing  that  can  be  done  with  the  ideas  here  is  the  following.  Suppose  that  one 
wanted  to  control  a  satellite  to  rotate  stably  about  its  middle  axis,  and  that  this  is  done  by  means  of 
the  techniques  of  this  paper,  say  through  internal  rotors.  Then  one  uses  geometric  phases  to 
reorient  the  body.  One  way  to  do  this  is  to  relax  the  control  so  the  body  becomes  unstable,  then  let 
it  go  unstable  so  it  will  swing  around  its  homoclinic  orbit,  and  when  it  comes  back  to  the  opposite 
saddle  point,  one  imposes  the  stabilizing  control  again.  One  can  also  help  initiate  motion  along  the 
homoclinic  orbit  using  a  linear  control  near  the  saddle  point,  as  in  Bloch  and  Marsden  [1989a]. 
Some  of  these  ideas  are  reminiscent  of  those  in  Beletsky  [1981],  p.  72.  This  maneuver  achieves  a 
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reorientation  of  180°  about  the  body’s  long  axis,  but  one  can  imagine  reorienting  it  using  more 
general  phases.  See  also  §5. 


§1  The  Rigid  Body  with  an  External  Torque 


The  rigid  body  equations  with  a  single  torque  about  the  minor  axis  are  given  by 

h~h 


“l  =  I 

I3  -  I, 

h  ~  l2 

0)3  =  J - ~  Ci)jC02  +  u 


(1.1) 


where  the  Ij  are  the  principal  moments  of  inertia  and  we  assume  for  the  moment  that  Ij  >  I2  >  I3. 
There  are  no  essential  changes  in  the  analysis  if  the  torque  is  taken  about  the  major  axis.  Now  let 
us  implement  the  feedback  torque 


I,-I2 

u  =  -  e-^j — -  0)jCd2  . 


(1-2) 


Making  the  transformation  to  the  classical  momentum  variables 

=  IjOOp  i  =  1,2,  3, 


the  equations  of  motion  become 


where 


Io-I 


m]  =  aj  rr^mj 


m3  =  a3(l-e)m1m2 


1,-1, 


(1.3) 


1,-1, 


a _2 tl  o  _  ii  and  a  1 — 

al  “  I2I3  ’  1^3  *  “a  ®3  "  Ijl 


Remarkably,  this  feedback  system  (1.3)  has  two  constants  of  motion. 


Lemma  1.1  The  functions 
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and 


hf 


_irsL  +  ”L  +  "L_!_') 

2  v  Ij  I2  I3  1  -  e  / 


Mp  =  |(mj(l-e)  +  m|(l-e)  +  m|). 


are  constants  of  the  motion  ( we  assume  £  *  1  for  the  moment). 


(1.4) 

(1.5) 


The  proof  is  a  straightforward  computation.  The  form  of  the  constants  of  motion  is 
chosen  to  give  us  the  following 

Proposition  1.2  The  controlled  system  (1.3)  is  Lie-Poisson  with  Hamiltonian  HF  with 
respect  to  the  Lie-Poisson  structure 

{F,  G}F  =  -VMp-(VFxVG)  (1.6) 

This  follows  from  checking  that  m;  =  [nij,  Hp]p,  i  =  1,  2,  3  (For  a  discussion  of  Lie- 
Poisson  systems  see,  for  example,  Krishnaprasad  [1985]  and  references  therein). 

There  are  three  other  canonical  Poisson  structures  associated  with  the  system  (1.3),  which 
are  discussed  in  Bloch  and  Marsden  [1989b].  Note  also,  as  pointed  out  in  the  latter  paper,  that  for 
£  <  1,  the  invariance  group  of  the  Lie-Poisson  bracket  (1.6)  is  SO  (3),  while  for  e  >  1,  it  is 
SO  (2,1).  Thus,  we  have  a  deformation  of  the  structure  with  the  parameter  e.  (Deformations  of 
Poisson  structures  are  discussed  in  Weinstein  [1983]).  In  fact,  we  can  show  that  (1.3)  are 
generalized  rigid  body  equations.  The  generalized  rigid  body  is  discussed  in  Abraham  and 
Marsden  [1978]  and  in  Arnold  [1978];  for  e  >  1,  we  get  an  indefinite  rigid  body,  which  is 
discussed  in  the  work  of  Klein  [1897]. 

Theorem  1.3  The  equations  (13)  are  the  generalized  rigid  body  (Euler)  equations  for  the  Lie 
group  G£  =  {  A  |  A  €  SL(3) ,  ATZA  =  2),  where  1  is  the  quadratic  form  given  by  I  = 

diag((l  -  e)-1,  (1  -e)-1, 1). 


Proof  We  use  the  fact  (see,  for  example,  Marsden,  et.  al.  [1983])  that  any  Lie-Poisson  system  is 
the  reduction  of  the  system  on  T*G  obtained  by  declaring  the  Hamiltonian  to  be  left  invariant  on 
the  corresponding  group,  in  this  case,  G  =  G£.  Note  that  (1.6)  is  the  Lie-Poisson  bracket  (with  a 
minus  sign,  hence  we  take  a  left  invariant  extension)  for  the  group  G£  with  Mp  a  Casimir.  One 

can  check  explicitly  that  the  equations  may  be  written  in  Eulerian  form  as 
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where  Az  is  the  Poisson  tensor: 


m  =  AIVHF(m) 


0 


-m. 


m3  0 

-  e)m2  (1  -  E)mj 


(1  -  E)m2 

-(1  -  E)mj 
0 


Bloch  and  Marsden  [1989b]  analyzed  stabilizability  of  the  free  rigid  body  system  about  its 
intermediate  axis  using  the  feedback  u  =  -  Ea3m1m2  .  This  is,  of  course,  equivalent  to  analyzing 

stability  of  the  closed  loop  system  (1.3). 

Note  first  that  (1.3),  linearized  about  the  (relative)  equilibrium  (irq,  m2,  m3)  =  (0,  M,  0), 
has  one  zero,  one  stable  and  one  unstable  eigenvalue  for  e  <  1,  but  has  one  zero  and  two 
imaginary  eigenvalues  for  £  >  1.  The  Energy-Casimir  method  shows  that  the  system  is  (non- 
linearly)  stable  for  £  >  1.  For  e  =  1  the  system  is  gyroscopically  stable. 

Recall  that  the  Energy-Casimir  method  (see,  for  example,  Holm,  Marsden,  Ratiu  and 
Weinstein  [1985]  and  Krishnaprasad  and  Marsden  [1987])  involves  finding  a  constant  of  motion 
for  the  system  of  the  form  H  +  C,  where  H  is  usually  the  energy  and  C  is  a  Casimir  function 
or  a  Casimir  function  plus  a  momentum  map,  such  that  H  +  C  has  a  critical  point  at  the  relative 
equilibrium  of  interest.  Definiteness  of  the  second  variation  of  H  +  C  at  the  critical  point  is 
sufficient  to  prove  stability  (at  least  in  finite  dimensions — the  infinite  dimensional  case  is 
discussed  in  the  cited  references).  In  the  case  here,  we  used  the  energy  Hp  and  the  Casimir  Mp  . 
It  is  easy  to  check  that  Mp  is  indeed  a  Casimir — it  commutes  with  every  dynamical  variable  under 
the  Poisson  bracket  { ,  }p.  We  omit  the  proof  here,  but  state  the  theorem. 


Theorem  1.4  The  rigid  body  equations  with  a  single  torque  about  the  minor  axis  and  with 
feedback  u  =  -  ea3mj  m2,  i.e.,  the  system  (1.3),  is  stabilized  about  the  relative  equilibrium  (m]5 
m2,  m3)  =  (0,  M,  0)  for  e  >  1. 


§2  Recovery  of  the  Externally  Torqued  System  and 
the  Heavy  Rigid  Body 

As  discussed  in  the  introduction,  the  fact  that  the  system  of  §1  with  external  torque 
feedback  has  conserved  quantities  and  a  Lie-Poisson  structure,  leads  one  to  ask  if  there  is  a 
mechanical  extension  of  the  system  to  larger  system  where  the  closed  loop  dynamics  is  realized  by 
an  internal  torque  feedback.  In  this  section  we  show  that  such  a  mechanical  extension  does  indeed 
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exist — a  rigid  body  carrying  three  symmetric  rotors  with  associated  internal  torques.  The  Lie- 
Poisson  structure  and  stability  of  this  system  was  first  analyzed  by  Krishnaprasad  [1985]  and  also 
by  Sanchez  de  Alvarez  [1986]. 

Consider  a  rigid  body  carrying  one,  two  or  three  symmetric  rotors.  Denote  the  system 
center  of  mass  by  0  in  the  body  frame  and  at  0  place  a  set  of  (orthonormal)  body  axes.  Assume 
that  the  rotor  axes  are  aligned  with  principal  axes  of  the  carrier  body. 

The  configuration  space  of  the  system  is  SO  (3)  x  S1  x  S1  x  S1  with  tangent  space  denoted 
T(S0(3)  x  S1  x  S1  x  S1). 

Let  be  the  inertia  tensor  of  the  main  body,  Irolor  the  diagonal  matrix  of  rotor 

inertias  about  the  principal  axes  and  1^.'  the  remaining  rotor  inertias  about  the  other  axes.  Let 
Ijock  =  I  body  +  I,.olor  +  I^tor  be  the  locked  inertia  tensor  (i.e.,  with  rotors  locked)  of  the  full  system; 
this  definition  coincides  with  the  usage  in  Marsden  et.  al.  [1989]. 

The  Lagrangian  (kinetic  energy)  of  the  free  system  is  the  total  kinetic  energy  of  the  body 
plus  the  total  kinetic  energy  of  the  rotor,  i.e., 

L  =  |  mbo<Iyn)  +  1  £2  £2  +  \  (£2  +  £2,)-  Iro„r(£2  +  £2,)  (2.1) 

=  5  (“  Olock  -  WJO)  +  5  (O  +  £V'W<o  +  o,>  <2.2) 

where  £2  is  the  vector  of  body  angular  velocities  and  £2r  is  the  vector  of  rotor  angular  velocities 
about  the  principal  axes  with  respect  to  a  body  fixed  frame. 

By  the  Legendre  transform,  the  conjugate  momenta  are: 

3L 

»  -  ^  -  <W  -  W  o  +  w<0  +  o,)  -  W  n  +  wo,  (2.3) 
W»<o  +  at)  (2.4) 

and  tlie  equations  of  motion  with  internal  torques  (controls)  u  in  the  rotors  are 

m  =  m  x  Q  =  m  x  (Ilock  -  I^r^m  -  /)  (2.5) 

/  =  u.  (2.6) 

We  now  show'  how  to  recover  the  feedback  equations  with  a  single  external  torque  from 
§1.  To  do  so,  it  is  convenient  to  use  the  variables  7t  =  m  -/  and  /.  Let  I  =  Ilock-  Srolor  (for 
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explicit  computation  we  will  let  I^y  =  diag(Ij,  I2, 13)  as  before  and  let  I  -  Ilock  ^otor  “ 
diagOj,  I2,  T3)).  Then  we  get  the  equations 


7t  =  (7t  +  /)  X  I  *71  -  U 

(2.7) 

l  =  U. 

(2.8) 

Then  we  have 

Theorem  2.1  There  is  a  choice  of  internal  torque  feedback  u(n,  l)  such  that  the  body 
dynamics  in  the  system  (2.7)-(2.8)  are  precisely  those  of  system  (1.3)  (with  external  torque 
feedback ) 


Proof:  Firstly  set 


u(tc,  1)  =  lx  r’rc - u'(7t), 


(2.9) 


so  the  system  (2.7)-(2.8)  becomes: 


Jt  =  JC  X  I  !7C  +  u'(7t) 


/  =  /  X  I  —  u'(7t)  . 


(2.10) 

(2.11) 


Now  suppose  that  I  =  diag(T j ,  T2, 13),  Tj  >  I2  >  T3  and 

u'(7t)  = 


0 

0 


L  ^2  J 


(2.12) 


Here  al5  a^  and  a3  are  defined  as  in  §1,  7t  =  (Ttj,  k2,  tc3)t,  and  l-(lvl2,l  3)T-  Thten  the 
equations  (2.7),  (2.8)  reduce  to 

TCj  — 

=  **2^3^  1 


—  a^(l  ^2 
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•  _  Iff  3  2 

*i  ~  ~ 

I3  I2 

i  _  _  ^3 

*2  “  ~  ~ 

II  I3 

/3  =  —  (/j -i- e7Cj)  —  —  (/2 -+- E7C2),  (2.13) 

T2  I, 

in  which  the  it  equations  are  precisely  (1.3).  ■ 

In  fact,  we  see  from  (2.10)  and  (2.1 1)  that  we  can  realize  any  external  torque  feedback  for 
the  rigid  body  with  our  rigid  body  plus  rotor  system. 

Consider  now  the  following  interesting  case.  Let  m  =  n  +  /  as  before  and  take  as  controls 
u'(ji,  m).  Then  (2.10)  and  (2.1 1)  may  be  written 

7t  =  7t  x  I_17t  +  u'(7t,  m) 

m  =  m  x  I_17t .  (2.14) 


In  particular  we  have 

Proposition  2.2  For  u'(7t,  m)  =  -  x  m  where  C  is  a  constant,  and  %  is  o.  constant 
(body  fixed)  vector,  the  equations  (2.14)  are  the  equations  for  the  heavy  top.  Hence  these  admit 
the  conserved  Hamiltonian  H=  1 7tT_17t  +  Cx  m  and  kinematic  conserved  quantities 

(casimirs),  Cj  =  ||m||2  and  =  Tt  m. 

This  may  be  verified  by  referring  to  the  heavy  top  equations  in  Marsden,  Ratiu  and  Weinstein 
[1984],  for  example.  The  equations  are  Lie-Poisson  on  the  dual  of  the  Euclidean  Lie  algebra  e(3). 

§3  The  Hamiltonian  Structure  of  the  Rigid  Body  with  Three  Rotors 

under  Feedback 

The  equations  of  motion  for  the  three  rotor  system  with  internal  torque  controls  are 


m  =  m  x  £2  =  m  x  (Ilock  -  Irotor)  !(m  -  7) 


(2.5) 
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i  =  u 

We  haive  the  following  result: 


(2.6) 


Theorem  3.1  For  the  feedback 

u  =  k(m  x  (Ilock  -  “ 0)«  (3- 1 ) 

where  k  is  a  constant  real  matrix  such  that  the  matrix  J  =  (1-  k)-1(Ilock  - 1^^)  is  symmetric, 
the  system  (2.5)  reduces  to  a  Hamiltonian  system  on  so (3)*  with  respect  to  the  standard  Lie- 
Poisson  structure  {F,  G}(m)  =  -  m  (VF  x  VG) 

Proof  We  have 

l  =  u  =  km  (3.2) 

which  is  an  acceleration  feedback.  (Note  that  also  u  =  k((I)ockQ  +  Irolor^T)  x  Q).  Therefore, 

km  -  /  =  p,  (3.3) 

where  p  is  an  arbitrary  constant.  Hence  the  closed  loop  system  becomes 

m  =  m  x  (Ilock  -  Irotor)“1(m  -  /) 

=  m  x  diock  -  -  km  +  p) 

=  mxa^-Wkl-kXm-^)  (3.4) 

where  £  =  -  (1  -  k)-1p.  Under  the  hypothesis  on  k,  we  can  define  the  k-dependent  "inertia 
tensor" 

i  -  (i  -  kr'fltod  -  w>  <3-5> 


Then  the  equations  become 


where 


m  =  VC  x  VH 
C  =  \  ||m||2 

h  =  \  (m  -  -  %) . 


(3.6) 

(3.7) 


and 


(3.8) 
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Equations  (3.6)  are  Hamiltonian  on  so(3)*  with  respect  to  the  usual  Lie-Poisson  structure. 


Remark  1  Note  that  (3.3)  is  equivalent  to 

k^lock  ^  +  4otor  ~  ^rotor  ^  +  =  P 

i.e., 


(^rotor  ^  ^rotor  )^r  “  ^^lock  ^rotor)  ^  P 


Thus  for 


Hock  "  ^rotor 


(3.9) 


one  specializes  to  the  usual  dual  spin  case  (see  Krishnaprasad  (1985)  and  Sanchez  de  Alvarez 
(1986))  where  the  acceleration  feedback  is  such  that  each  rotor  rotates  at  constant  angular  velocity 
relative  to  the  carrier.  Also  note  that  the  Hamiltonian  in  (3.8)  is  indefinite  for  high  gains. 


Remark  2  Furthermore,  if  we  set  mb  =  m  -  the  equations  become 

mb=  (mb+  £)  x  I-1  mb 


(3.10) 


which  is  a  k-dependent  form  of  the  dual-spin  equations.  Hence  one  can  apply,  for  example,  the 
Morse-theoretic  analysis  of  Krishnaprasad  and  Berenstein  [1984]  to  this  situation.  For  a  related 
stability  analysis  using  the  energy-momentum  method  (see  L.S.  Wang  [1990]).  ♦ 


It  is  instructive  to  consider  the  case  where  k  =  diag(kj,  k2,  k3).  Let  I  =  (Ilock  -  ^otor)  = 
diag(fj,  I2,  I3)  and  the  matrix  I  satisfies  the  symmetry  hypothesis  of  Theorem  3.1 .  Then  /;  = 
p;  +  kjtnj,  i  =  1,  2,  3,  and  the  equations  become  m  =  m  x  VH  where 


l  f  ((1  -  k^m,  +  pj)2 
H  (1  -  kj)Tj 


(d  -  k2)m2  +  p2)2  ((1  -  k3)m3  +  p3)2  N 

(1  -  k2)I2  (1  -  k3)T3 


(3-11) 


It  is  possible  to  have  more  complex  feedback  mechanisms  where  the  system  still  reduces  to  a 
Hamiltonian  system  on  so(3)*.  If  we  set  u  =  k(m,  l)m,  the  key  to  reduction  is  that  k(m,  /)  m  - 
/  =  0  is  integrable.  The  natural  case  to  consider  is  k  =  k(m).  Then  we  have 
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Theorem  3.2  The  system  (2.5),  (2.6)  with  u  =  k(m)m,  where  k(m) 


kjM1 

k2(m)T 

k3(m)T 


and 


kj(m)  =  Vfyfm)  for  some  smooth  fy,  and  such  that  (Ilock  -  Irotor )  ^(m)  is  symmetric,  is 
Hamiltonian. 


Proof  We  have 

ki(rn),ih  -  /j  =  0  =  V0i(m)m -f  , 

i.e.,  5  [<t>i(m)  -  /J  =  0 

or  4>i(m)  -  h  =  Pi  • 

3<j) 

Thus,  k(m)  =  — ,  where  <3>  is  the  vector  potential  [<J>,  <{>2  <j>3]T  gives  integrability.  Then  (2.5) 
dm  ’  ’ 

and  (2.6)  reduce  to 

m  =  mx  (Ilock  - 1,0^  )_1(m  -  4>(m)  +  p) .  (3.12) 

This  is  Lie-Poisson  on  50(3)*  if  and  only  if  there  exists  an  H(m)  such  that 

VH  =  Ollock  “  ^rotor  “  *(“)  +  P)- 

By  equality  of  mixed  partial  derivatives,  this  holds  if 

“  (erdiock  "  W  "  ®(m)  +  p))  =  ~  (ej  dlock  -  W  )_1(m  -  O(m)  +  p)), 

where  e;  is  the  i*  unit  vector.  This  is  true  if  and  only  if  e;-(I]ock  -  Irotor  )-1k(m)ej  =  ej 
aloCk  -  W  )_1k(m)  Cj ,  i.e.,  (Ilock  -  \otor  r1k(m)  is  symmetric.  ■ 

To  compute  H,  we  need  to  solve  the  set  of  partial  differential  equations 
VH  =  (Iloek-  Irotor  )-1(m  -  O(m)  +  p). 

Consider  the  special  case  where  (Ilock  -  Irolor )  and  k(m)  are  diagonal,  (Ilock  - ) 

„  _  _  B0> 

=  diag(Ij,  I2,  I3)  and  k(m)  =  diag  (k^,  k^,  k3m3).  Then  k(m)  =  —  (m)  where 
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Hence 


<E>(m) 


1 

2 


klm? 

k2m2 

k3m3 


0H 

dmj 


k;m2 

_ i 


+  Pi 


i  =  1,  2,  3. 


3 

Thus,  letting  H  =  ^  H  ^ro)  we  have 

i=  l 


h 


(  m2 


k;m: 

j _ i 


3  \ 

+  Pimi 


(3.13) 


This  procedure  clearly  generalizes  to  k^ir^)  an  arbitrary  polynomial  in  nr,  yielding  a 
large  class  of  Hamiltonian  feedback  systems. 


§4  The  Rigid  Body  with  a  Single  Rotor- 
Feedback  Stabilization  by  the  Energy-Casimir  Method 


We  now  consider  the  equations  for  a  rigid  body  with  a  single  rotor.  This  reveals  the  essentials  of 
the  dynamics  and  Hamiltonian  structure  and,  further,  we  are  able  to  obtain  the  result  that  with  a 
single  rotor  about  the  third  principal  axis,  quadratic  feedback  with  sufficiently  high  gain  stabilizes 
the  system  about  its  intermediate  axis  -  a  result  similar  to  that  of  the  single  gas  jet  (external  torque) 
case. 

We  adopt  the  following  notation:  Let  the  rigid  body  have  moments  of  inertia  Ij  >  ^  >  I3  as 
before  and  suppose  the  symmetric  rotor  is  aligned  with  the  third  principal  axis  and  has  moments  of 
inertia  Jj  =  J2  and  J3.  Let  a>j,  i  =  1,  2,  3,  denote  the  carrier  body  angular  velocities  and  let  a 
denote  that  of  the  rotor  (relative  to  a  frame  fixed  on  the  carrier  body).  Let 

diag  (Xlf  ^2,  X3)  =  diag  (J}  +  Ip  J2  +  J3  +  I3)  (4.1) 

be  the  locked  inertia  tensor.  Then  from  (2.3)  and  (2.4),  the  natural  momenta  are 


mi  =  (Ji +  »  i  =  1.  2 
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m3  “  ^3®3  +  h 


/3  =  J3(co3  +  a) 


(4.2) 


From  (2.5)  and  (2.6),  the  equations  of  motion  are: 


mx  = 


1  _  1 
\3  ^2  J 


m2  —  mim3 


m3  =  m1m2 


1  1 


f 


hmi 


1  1 


VA2 


/3  =  u 


(4.3) 


Choose  u  =  ka3m]m2  where  a3  = 


(k2  ^1/ 


.  Then  we  have 


Theorem  4.1  With  u  =  ka1m1m2  the  equations  (4.3)  reduce  to  the  equations 


rrij  =  m2 


^(1  -  k)m3  -  p  ^ 


= 


-m, 


(1  -  k)m3  -  p  I  m3m1 


m3  =  a3m1m2 


(4.4) 


which  are  Hamiltonian  on  so (3)*  with  respect  to  the  usual  Lie-Poisson  structure  with 
Hamiltonian 


^+^|  +  «1-k)m3-p)2' 
^2 


(1  -  k)I3  J 


J30  ~  k) 


(4.5) 


where  p  is  a  constant 
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Proof  First  note  that  p  =  /3  -  km3  is  conserved.  Eliminating  l3  gives  rise  to  equations  (4.3) 
which  are  Hamiltonian  with  respect  to  the  standard  bracket  {f,  g)(u)  =  -  m  (Vf  x  Vg)  and  the 
given  H.  ■ 


Remark  When  k  =  0,  we  get  the  equations  for  the  free  rotor.  We  get  the  dual  spin  case  J3a  = 

0,  when  k  Ilock  =  Irotor  from  (3.9),  or,  in  this  case  k  =  —  .  One  can  check  that  for  this  k,  p  = 

X3 


(l-k)a,  a  multiple  of  a.  ♦ 


We  can  now  use  the  Energy-Casimir  method  to  prove 

I, 

Theorem  4.2  For  k  >  1  —  —  (and  p  =  0)  the  system  (4.4)  is  stabilized  about  the  middle 

\2 

axis,  i.e.,  about  the  relative  equilibrium  (0,  M,  0). 

Proof  Consider  the  Energy-Casimir  function  H  +  C  where  C  =  cp(m2),  and  m2  =  m2  +  m2  + 
m2 .  The  first  variation  is 


S(H  +  C)  =  —  -  —  +  — -  2  +  - 1 — - — -  5m3  -f  9,(m2)(m18m1  +  m28m2  +  m38m3).(4.6) 

Xj  X2  a3 


This  is  zero  if 


m 


+  tp'nij  =  0, 


m? 

7  +  <Pm2  =  0, 
K1 


(1  -  k)m3  -  p 

j  +  <p  m3  - 


0 


(4.7) 


Then  we  compute 


(5m t)2  (Sn^)2  (1  -  k)(5m3)2 


52(H  +  C)  =  - — —  +  - — —  +• 
X  j  X2 


+  (p,(m2)((8m1)2  +  (Smj^-f  (8m3)2) 

+  9"(m2)(m18m1  +  n^Smj  +  m38m3)2  . 


(4.8) 
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If  p  =  0,  i.e.,  /3  =  km3,  (0,  M,  0)  is  a  relative  equilibrium  and  (4.7)  are  satisfied  if  tp' 
In  that  case, 


X2 


S2(H  +  C)  =  (8m  j)2 


_1__  J_ 
X2  j 


+  (8m3): 


1  -  k  1 


+  <p"(5m2)2 


1  1 

Now - - 2 


^  <  0  for  Ij  >I2  >I3  .  Clearly  for  k  sufficiently  large  — <  0, 

’  X0 


X,  X, 


/v2  /v  1  ^2  J  /V2 
so  if  we  choose  cp"  <  0  we  get  negative  definiteness  and  hence  stability.  (Note  that  for  k  =  0  the 

second  variation  is  indefinite  as  it  should  be.)  ■ 


Note  that  we  don't  get  precisely  the  gas  jet  equations  as  we  do  in  §2,  but  we  achieve  the 
same  effect  -  for  k  greater  than  a  certain  value  we  achieve  stability  about  the  middle  axis.  Also, 
we  are  using  a  much  simpler  feedback  than  in  Theorem  2.1  of  §2,  as  we  don't  need  to  "cancel" 
the  /3  dynamics.  Note  also  that  while  the  Hamiltonian  here  is  indefinite  for  k  sufficiently  large, 

the  conserved  momentum  here  is  the  usual  definite  momentum,  in  contrast  to  the  possibly 
indefinite  "momentum"  of  the  gas  jet  case  in  §1. 

Corresponding  to  the  Hamiltonian  (4.5)  there  is  a  (nonphysical)  Lagrangian  found  using 
the  inverse  Legendre  transformation 


(1  -  k)m3  -  p 

"  - h  - 

s  =  _  _ 2 _  (< 

(1  -  k)I3  (1  -  k)I3 

Note  that  5j,  o>2  and  6^  are  the  angular  velocities  0)j,0)2,  and  CtXj  for  the  free  system,  but 
that  S  is  not  equal  to  a.  In  fact, 


fi 


a  km3 

(1  -  k)  ~(1  —  k)J3  ’ 


(4.10) 
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Thus  the  equations  on  TS0(3)  determined  by  (4.9)  are  the  Euler- Lagrange  equations  for  a 
Lagrangian  quadratic  in  the  velocities,  so  the  equations  can  be  regarded  as  geodesic  equations.  The 
torques  can  be  thought  of  as  residing  in  the  velocity  shift  (4.10).  Using  the  free  Lagrangian,  the 
torques  appear  as  generalized  forces  on  the  right  hand  side  of  the  Euler-Lagrange  equations.  Thus, 
the  d'Alembert  principle  can  be  used  to  describe  the  Euler-Lagrange  equations  with  the  generalized 
forces.  However,  this  latter  approach  hides  the  useful  fact  that  the  equations  are  actually  derivable 
from  a  Lagrangian  (and  hence  a  Hamiltonian)  in  velocity  shifted  variables.  Also,  for  problems  like 
the  driven  rotor,  one  might  think  that  this  is  a  velocity  constraint  and  should  be  treated  by  using 
constraint  theory.  Again  for  the  particular  problem  at  hand,  this  can  be  circumvented  and  standard 
methods  are  in  fact  applicable,  as  we  have  demonstrated. 

§5  Phase  Shifts 

In  this  final  section  we  discuss  an  attitude  drift  that  can  occur  in  the  system  and  a  method 
for  correcting  it 

If  the  system  (3.10)  is  perturbed  from  a  stable  equilibrium,  and  the  perturbation  is  not  too 
large,  the  closed  loop  system  executes  a  periodic  motion  on  a  level  surface  (momentum  sphere)  of 
the  Casimir  function  ||  m^  +  £  ||2  in  the  body-rotor  feedback  system.  This  leads  to  an  attitude  drift 

which  can  be  thought  of  as  rotation  about  the  (constant)  spatial  angular  momentum  vector. 

We  will  calculate  the  amount  of  this  rotation.  This  can  be  done  following  a  method 
developed  by  Montgomery  [1990]  for  calculating  the  phase  shift  for  the  single  rigid  body.  Other 
pertinent  work  on  phases  may  be  found  in  Marsden,  Montgomery  and  Ratiu  [1990]  and 
Krishnaprasad  [1990]. 

Recall  that  the  equations  of  motion  for  the  rigid  body-rotor  system  with  feedback  law  (3.1) 
may  be  written  (see  3. 10)  as 

mb  =  (mb  +  %)  x  J-1mb 

=  VC  x  VH  (5.1) 


where  mb=  m  +  (1-  k)  !(km  - 1)  and  £  =  -  (1-  k)  !(km  - 1)  and 

C  =  2  Hmb+£||2  (5.2) 

H  *  \  mb'J'lmb  • 


(5.3) 
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Here  H  is  the  conserved  Hamiltonian  and  C  =  constant  determines  a  momentum  sphere  in  the 
reduced  phase  space  centered  at  t,  e  R3  =  so(3)*. 

The  attitude  equation  for  the  rigid-body  rotor  system  is 

A  =  AD  , 

where,  in  the  presence  of  the  feedback  law  (3.1), 

=  Oiock-WV-O 

=  ^lock- W)_1(m~km  +  P) 

=  J_1(m  -  4)  =  J-1mb  . 

Therefore  the  attitude  equation  may  be  written 

A  =  A(l_1mb)A  (5.4) 

where  A  denotes  the  canonical  map  from  vectors  in  R3  to  elements  of  so  (3), 

-Vl-|A  [“  0  -v3  v2 

v2  =  v3  0  -Vj 

’  V3  L-V2  Vj  0 

and  A  is  the  attitude  matrix,  A  e  S0(3). 

The  net  spatial  (constant)  angular  momentum  vector  is  now  given  by 

p  =  A(mbH)  (5-5) 

Then  we  have  the  following 

Theorem  5.1  Suppose  that  the  solution  of 

mb  =  (mb  +  4)  x  J-1mb  (5.6) 

sweeps  out  a  periodic  orbit  of  period  T  on  the  momentum  sphere,  |j  0^+  4 II2  =  il  (i  ||2,  enclosing  a 

solid  angle  $wlid.  Let  £2av  denote  the  average  value  of  the  body  angular  velocity  over  this 
period.  Let  E  denote  the  constant  value  of  the  Hamiltonian  and  let  ||  (i  ||  denote  the  magnitude  of 
the  angular  momentum  vector.  Then  the  body  undergoes  a  net  rotation  A0  about  the  spatial 
angular  momentum  vector  (I  given  by  the  formula 
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A0 


2ET  (  T 
II 4 II  +  Hull 


+  ^solid 


(5.7) 


Proof  Consider  the  reduced  phase  space  (the  momentum  sphere) 

P„  =  {mj  ||mb+^p  =  |Ui|pj,  (5.8) 

for  (J.  fixed.  Suppose  we  have  a  periodic  orbit  in  of  period  T.  We  then  compute  how  A(t) 
changes.  Now 

mb(to  +  T)  =  m^)  (5.9) 

and  from  momentum  conservation 

H  =  A(T  +  t0)(mb(T+I0)  +  y 

=  AdoKmtdo)  +  Q  .  (5.10) 

Hence 

A(T  +  t0)A(t0r 1  it  =  p . 


Thus  A(T  +  t0)A(t0)  1  is  an  element  of  G^,  the  subgroup  of  the  group  of  rotations  that  fixes  p, 
so 


A(T  +  t0)A(t0)_1 


exp 


A0 


\ 


_P_ 

Mj 


(5.11) 


for  some  A0,  which  we  wish  to  determine. 

Consider  a  phase  trajectory  of  our  system 

z(t)  =  (A(t),  mb(t)),  z(tQ)  =  z0  (5.12) 

with  A(to)  =  I  and  m^tg)  =  p  -  £;,  the  body  thus  being  in  the  reference  configuration  at  t  =  tg. 
The  two  curves  in  phase  space 

Cx  =  (z(t) :  Iq  <  t  <  tQ+  T}  (the  dynamical  evolution  from  Zg) 


and 


r 

(  > 

jexp 

0  — 

Hull  J 

zo:O<0<A0f  intersect  at  t  =  T. 


Thus  C  =  Cj  -  C2  is  a  closed  curve  in  phase  space  and  from  Stokes'  theorem  we  have, 
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Jc  P<iq  "Jc  Pdq  =  d(pdq)  (5.13) 

3 

where  pdq  =  ^  Pjdqj  and  where  q;  and  pj  are  configuration  space  variables  and  conjugate 
i=  1 

momenta  in  the  phase  space  and  X  is  a  surface  enclosed  by  the  curve  C. 


Evaluating  each  of  these  integrals  will  give  the  formula  for  A8.  To  see  this,  first  consider 
the  dynamical  evolution  along  Cj.  Letting  co  be  the  spatial  angular  velocity,  we  get 

=  n-o)  =  A(JQ  +  £)-Aft  =  (K2 +  £)•£} 


=  JQ Q  +  £ Q  . 

Hence 

ic,pdq  =  Jc,  p' *  dl  ‘  Jo  ,Q£2dt  +  Jo  ^dt 


=  2ET+ft-n„)T 

since  the  Hamiltonian  is  conserved  along  orbits.  Similarly,  along  C2  we  have, 


where  we  substituted  the  spatial  angular  velocity  along  C^. 


(5.14) 


(5.15) 


(5.16) 


Finally  we  note  that  the  map  7t^  from  the  set  of  points  in  phase  space  with  angular  momentum  (I 
to  satisfies 

Jjz  d(pdq)  -  fjn  (Z)dA  =  ±  Hull  «>«*,  (5.17) 


where  dA  is  the  area  form  on  the  two-sphere  and  rc^(Z)  is  the  spherical  cap  bounded  by  the 
periodic  orbit 


{mb(t)  I  to^t^to+T}  cP^. 
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From  (5.15),  (5.16)  and  (5.17)  we  have  the  result.  ■ 

Remark  5.2  a  When  £  =  0,  (5.7)  reduces  to  Montgomery's  formula  [1990], 

b  This  theorem  may  be  viewed  as  a  special  case  of  a  scenario  that  is  useful  for  other  systems, 
such  as  rigid  bodies  with  flexible  appendages.  Phases  are  usually  viewed  as  occurring  in  the 
reconstruction  process,  which  lifts  the  dynamics  from  to  J_1(p),  where  J  :  P  ->  g*  is  the 
momentum  map  for  a  mechanical  system  with  symmetry  and  P..  =  J_1(p)/Gu  is  the  reduced  space; 
see  Marsden,  Montgomery  and  Ratiu  [1990].  The  cotangent  bundle  reduction  theorem  states  that 
P^  itself  is  a  bundle  over  T*S,  where  S  =  Q/G  is  shape  space.  The  fiber  of  this  bundle  is  O^, 
the  coadjoint  orbit  through  p.  For  a  rigid  body  with  three  internal  rotors,  S  is  the  three  torus  T3 
parametrized  by  the  rotor  angles.  Controlling  them  by  a  feedback  or  other  control  and  using  other 

conserved  quantities  associated  with  the  rotors  as  we  have  done,  leaves  one  with  dynamics  on  the 
"rigid  variables"  O^,  the  momentum  sphere  in  our  case.  Then  the  problem  reduces  to  that  of 

lifting  the  dynamics  on  to  J-1(p)  with  the  T*S  dynamics  given.  For  G  =  S0(3)  this 

"reduces"  the  problem  to  that  for  geometric  phases  for  the  rigid  body  given  by  Mongomery  [1990]. 

♦ 

Finally,  following  Krishnaprasad  [1990]  we  show  that  in  the  zero  total  angular  momentum  case 
one  can  compensate  for  this  drift  using  two  rotors. 

The  total  spatial  angular  momentum  if  one  has  only  two  rotors  is  of  the  form 

p  =  A(IlockQ  +  bjdj  +  b2a2)  (5.18) 

where  the  scalars  dj  and  d2  represent  the  rotor  velocities  relative  to  the  body  frame.  The 
attitude  matrix  A  satisfies 

A  =  Afl,  (5.19) 

as  above.  If  p  =  0,  then  from  (5.18)  and  (5.19)  we  get, 

A  =  -AWii'k  b,)A a,  +  ( I,^k  bj)'  6 t2)  .  (5.20) 

It  is  well  known  (see  for  instance  Brockett  [1973]  or  Crouch  [1986])  that  if  we  treat  the  a.,  i  =  1, 
2  as  controls,  then  attitude  controllability  holds  iff 
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(U  and  (Ij^?kb2)A  generate  so(3),  (5.21a) 

or  equivalently,  iff  the  vectors 

^lock  bi  and  Ij^?k  b-,  are  linearly  independent  (5.21b) 

Moreover,  one  can  write  the  attitude  matrix  as  a  reverse  path-ordered  exponential 

(5.22) 

A  key  observation  is  that  the  right-hand  side  of  (5.22)  depends  only  on  the  path  traversed  in  the 
space  T2  of  rotor  angles  (ap  a2)  and  not  on  the  history  of  velocities  or.  This  can  be  easily 
checked  by  carrying  out  a  time  rescaling  t  h4  p(t).  Hence  the  formula  (5.22)  should  be  interpreted 
as  a  "geometric  phase".  Furthermore,  the  controllability  condition  can  be  interpreted  as  a  curvature 
condition  on  the  principal  connection  on  the  bundle  T2  x  SO (3)  — »  T2  defined  by  the  so  (3)- 
valued  differential  1-form, 

e(ap  a2)  =  -(( I^k  b/  dttj  +  (l,^k  b2)A  da2).  (5.23) 

For  further  details  on  this  geometric  picture  of  multibody  interaction  see  Krishnaprasad  [1990]  and 
Krishnaprasad  and  Yang  [1990]. 


A(t)  =  A(0)P  exp 


-  L  ( bi>A  «i«») + b2)A  «2  ««  do 


Conclusions 

Natural  mechanical  systems  subject  to  exogenous  forces  determined  by  suitable  classes  of 
feedback  laws  can  sometimes  be  modeled  as  hamiltonian  systems.  Often  the  hamiltonian  structures 
so  derived  may  be  viewed  as  deformations  (by  feedback  gains)  of  the  hamiltonian  structure 
governing  the  open  loop  unforced  system.  The  methods  of  geometric  mechanics  such  as  reduction, 
reconstruction  phases,  energy-Casimir  and  energy-momentum  algorithms  for  stability  analysis, 
prove  to  be  naturally  applicable  to  such  feedback  systems.  The  present  paper  illustrates  this  via  a 
careful  study  of  the  problem  of  rigid  body  control  using  external  torques  (as  implemented  by  gas 
jets)  and  internal  torques  (via  reaction  wheels/rotors),  and  the  relationships  between  these  two 
methods  of  control.  We  reveal  a  Hamiltonian  structure  for  the  controlled  system  that  involves  a 
velocity  shift  (or  a  gyroscopic  term),  a  technique  that  should  be  of  general  utility.  We  especially 
note  that  we  find  a  nonobvious  but  rich  Hamiltonian  structure  for  the  controlled  system,  despite  the 
presence  of  controlling  torques. 
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One  of  the  main  accomplishments  of  the  present  work  is  the  stabilization  analysis  of  these 
examples. as  the  feedback  parameter  is  increased.  We  also  calculate  the  geometric  phase  (the 
attitude  drift)  in  the  presence  of  this  control  when  one  is  executing  a  periodic  motion  in  the  carrier 
body  angular  momentum  phase  space. 

A  general  theory  for  feedback  systems  of  gyroscopic  type  appears  in  the  recent  thesis  of 
Wang.  We  note  however  that,  in  the  examples  of  the  present  paper,  we  go  on  to  consider  feedback 
laws  that  also  effectively  alter  the  metric  structure  underlying  a  simple  mechanical  system.  The 
introduction  of  geometric  phases  (a  subject  of  great  current  interest  to  physicists)  in  the  analysis  of 
attitude  drift  due  to  perturbations  from  relative  equilibria  is  one  of  the  novel  features  of  the  paper. 
The  methods  we  use  for  both  the  phase  and  stability  calculations  should  be  applicable  to  more 
complex  many-body  systems  such  as  those  incorporating  flexible  elements. 

In  work  under  way,  we  are  also  considering  problems  of  instability  in  feedback  systems.  It 
appears  to  be  possible  to  establish  instability  theorems  based  on  second  variation  criteria  analogous 
to  the  classic  work  of  Kelvin  and  Chetayev.  Details  will  appear  in  a  forthcoming  paper. 
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